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Abstract. Phase operators and phase states are introduced for irreducible representations of 
the Lie algebra 5u(3) using a polar decomposition of ladder operators. In contradistinction with 
su(2), it is found that the 5u(3) polar decomposition does not uniquely determine a Hermitian 
phase operator. We describe two possible ways of proceeding: one based in imposing SU (2) 
invariance and the other based on the idea of complementarity. The generalization of these 
results to SU (n) is sketched. 



1. Introduction 

Phase is a unique concept for the proper understanding of classical optical phenomena. It 
is therefore surprising that, at a foundational level, quantum optics can apparently subsist 
without a quantum phase. One can use, for example, the better behaved field-quadrature 
operators 1 1 1, or be content with a pragmatic approach in which phase is a parameter that can 
be efficiently estimated |2-7|. It is equally possible to represent states as quasidistribution 
functions in phase space and specify their phase properties by classical angles [8 - 13J. Finally, 
it is also entirely reasonable to approach the problem from an operational perspective (14 -18 1 



that emphasizes the apparatus involved in measurements on the system; the phase then refers 
to a feature of this apparatus. 

However, if one adheres to the orthodox interpretation of quantum mechanics and regards 
phase as a physical property then surely it ought to be represented by a Hermitian operator. In 
other words, the phase variable should be subject to quantization and, for a sufficiently small 
number of particles, quantized phase effects should be accessible to experimentation. 

Despite the difficulties borne out of the first and eminent attempts p9}j2l1 |, and ultimately 
ascribed to the semiboundedness of the eigenvalue spectrum of the number operator, 
significant progress has been achieved in the last fews years in clarifying the status of a 
quantum phase operator. The primary objective in this context has been the description of 
the phase of a single-mode field, or, equivalently, of a harmonic oscillator. The progress made 
is manifest and the work on this subject has already been reviewed |22}j25| . 

Although the definition of the absolute phase is in itself an interesting problem, such 
an absolute phase has no meaning from a practical point of view. Strictly speaking, only 
relative variables are of interest in physics. Most, if not all, methods of phase measurement 
are arrangements determining the relative phase between two different modes. One might 
falsely expect that the relative phase should be constructed merely as the difference of phases. 
Perhaps surprisingly, experience demonstrates that this is not the case: there are theoretical 
and experimental results that cannot be accounted for using the difference of phases (26] . In 
view of this, one should start a study of the relative phase without any previous assumption 
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about single-mode phases. In particular, the conjugate variable to a relative phase is a number 
difference that is not bounded from below. Thus, it is reasonable to expect that the relative 
phase will be free of the problems arising in the one-mode case. 

In the characterization of the relative phase for two-mode fields, the Stokes parameters 
play an important role | [27| . From an experimental point of view, they are measurable 
quantities. An essential observation is that they are formally also elements of 5u(2), which 
turns out to be the dynamical symmetry algebra of a qubit, in the modern parlance of quantum 
information |28|. Interesting links with finite quantum systems have been discussed in 
references [29, 30|. In fact, as shown time ago in the pioneering works of Levy-Leblond |3T| 
and Vourdas |32-34| (see also |35 1), the polar decomposition of 5u(2) operators gives rise 
to a bona fide phase operator, which is also complementary to the population difference. 



generalizing somehow Dirac original idea p6}j38| 



One might believe — again falsely — that the passage from two-level systems and 5u(2) 
to three-level systems and 5u(3) would be immediate. We realize in this paper that this is not 
so: there appears to be no general phase operator for 5u(3) that simultaneously verify a polar 
decomposition, hermiticity and adequate commutation relations. 

The polar decomposition of an operator is always possible in any dimension (although 
it is generally not unique). On the other hand, complementarity in finite dimensional spaces 
is usually implemented via finite Fourier transformations. It will be shown that the happy 
coincidence where both concepts occur in the same problem must be in general abandoned 
for 5u(3) and more generally forsu(^): barring exceptional circumstances, complementarity 
and polar decompositions are apparently incompatible. 

Because it is generalized much more easily, the definition of phase operators obtained via 
polar decompositions will be studied in details, with special attention to the Lie algebra 5u(3). 
For the general 5u{n) case, it is natural to define n—l relative phase operators; we find that 
they do not in general commute, except in very specific circumstances where the dimension of 
the system is n or goes to infinity. It will be shown how to quantify this lack of commutativity 
and how this can be related to simple counting arguments based on the geometry of 5u{n) 
weight space. 

We have found that considerable insight in the structure of phase operators, emphasizing 
the connection between polar decomposition methods, the abstract operators and the 
geometrical nature of the relative phase, is gained by introducing a coherent- state realization 
of the generators (see references [ (39l[40| for variations on this theme). These realizations also 
provide very useful calculational simplifications, particularly as representations become large 
and as we increase the rank ^ — 1 of 5u{n). We introduce this representation first for 5u(2), 
in section [2j reserving a wealth of mathematical details for [Appendix A In particular, these 



realizations allow us to reach some conclusions about the commutativity of phase operators 
for 5u{n) with ^ > 3, in the limit of large representations. This limit can be seen as a classical 
limit, thus allowing our conclusions to be checked against classical concepts associated with 
relative phases. 

2. SU (2) phase operators 

The complex extension of the 5u(2) algebra is generated by the operators {h^e^^i-} with 
commutation relations 

[h,e±] = ±e±, [e^,e-]=2h. (2.1) 

In the (2 j + 1) -dimensional space Vj spanned by the vectors {\jm) : m = —7, . . . ,7}, which 
is the carrier of the irreducible representation (irrep) with spin 7, the generators act in the 
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standard way: 



h\jm) =m\jm), e±\jm) = ^J{j^m){j ±m^\)\jm±\) . (2.2) 
The state = \jj) is the highest weight of the irrep, so that 

h\Xj)=j\Xj). e^\Xj)=^- (2.3) 

An advantageous reahzation providing a Hnk with the two-mode relative phase, is given 
by the Schwinger reahzation of 5u(2) in terms of two bosonic fields di and ^2 | [4T|[42| : 

ej^^d\d2^ e-^d\d\^ h ^ ^{d\d\ — d\d2) ^ (2.4) 

which are the building blocks of the Stokes operators |27|. They act on the two-dimensional 
harmonic oscillator basis |^i , ^2) related to the angular momentum states \jm) by ^1 +^2 = 2j 
and n\—n2 = 2m. 

If, for a moment we interpret di and ^2 as classical field amplitudes, it is apparent from 



(2.4) that the relative phase between the fields is encoded in e±. Therefore, it seems natural 
enough to look for a polar decomposition of the ladder operators |[43j 



ED, 



(2.5) 



where, using (2.2), D = y e_e- is 3. semi-positive self-adjoint operator and £" is a unitary 
operator that can be interpreted as the exponential of a putative phase operator. 

The rank of e- is one less than the dimension of so E is not completely specified. 
We remove the ambiguity by using cyclic boundary conditions so that E is the generator 
of an Abelian cyclic group |38| . The eigenvalues of E are then the quantized phases. As 
E'^j^^ = i, these eigenvalues are just O)^, where CO = exp[2;r// (2 j + 1)] and k = — j, ...J. 
Taking m modulo 27 + 1, we get in this way 



J 

n=-j 



(2.6) 



or, equivalently. 



/ 
1 



1 \ 



(2.7) 



/ 



The eigenstates of E are related to the basis states \jm) by a finite Fourier transform and are 
thus complementary to the \jm) states. 

Observe that E is not an SU (2) matrix. It can be written formally as the exponential of 
a Hermitian phase operator 9, but <p is not in general an element of the 5u(2) algebra. These 
observations hold (even in dimension 2 because of the choice of phase in completing E), even 
though both have well-defined actions on the basis elements. 

As heralded in the introduction, we now introduce as a tool of particular convenience 
a coherent state realization F for 5u(2). F is not Hermitian, although it is equivalent 
by similarity transformation to the "standard" Hermitian realization of p.2|), as shown in 



Appendix A| 

The SU{2) coherent states are defined by | [44}j46| 



(2.8) 
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where and Ry represent rotation about the z and y axes, respectively. To any vector |^) we 
associate the function 

\W) ^ W^i(p) = iXjlRyWR.m'V) . (2.9) 

Note that (<p) p is precisely the Husimi Q-f unction for the corresponding (pure) state |^). 
We can use the arbitrary nature of |^) to define the action of X G 5u(2) on ^T}{(p) by: 

X\^) ^ [r{X)^U(p) = {Xj\RyWRd(p)X\'i') ■ (2.10) 

Straightforward manipulations immediately produce the expressions 

d 



i^r(h) = -i- . 

dip 



(2.11) 

d 



e^^r{e^) = -(tant^)^ie±^'^ y^'dcp 

It is convenient to think of the coherent state 1 (jo) as localized around the coordinates (t^, (p) 
on the Bloch sphere. The relative phase is linked to the azimuthal angle (p while the parameter 
t3- is inessential for our purposes; to simplify the expressions for F, we choose t3- so that 
tan = — 1 and obtain 

d(p 



d(p 



One easily verifies that X r(X) preserves the commutation relations (2.1 



(2.12) 



and is thus 



a realization of 5u(2). F acts naturally in the infinite-dimensional space spanned by the 
exponential functions {e'^^ : 2m G Z} and equipped with scalar product 

{f\g)= r{(p)g{(p)d(p. (2.13) 

Jo 

The subspace of states with |m| < j is invariant under the action of F. The (normalized) 
basis elements of this invariant subsp ace are mapped to exponential functions \jm) ^ 
exp(/m(p) / >/27r and the action of ( 2.12 ) is 



T{h)\jm) =m\jm) r{e±)\jm) = {j^m)\jm±\) . (2.14) 
Under this inner product 

{jm'\T{e^)\jm) + ^n{\Y\e_)\jm) , (2.15) 



thus, the realization F is not Hermitian. However, we show in Appendix A how, for fixed j, F 
is equivalent to the standard Hermitian representation given in equation (2.2). 

The considerable merit of the realization F is that it is particularly well-suited to analyze 
the polar decomposition: the unitary matrix E in Vj is immediately obtained from the action 
of e~^^, the "phase" part of F(^_): 

E^,^ = F(£)^,^ = (jW|e-^'^|jm) , -j<m<j mod(2j + 1) .(2.16) 

Thus, simply shifts the basis state on the circle to its immediate neighbour e^^'""^)^, 
modulo 2 J + 1 . In addition, matrix elements of F becomes indistinguishable from those of the 
standard hermitian representation in the limit where m/ j ^ 0: this makes F also very well 
suited to analyze some limits of large representations, and analyze a transition between the 
quantum and classical phase. 
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In the basis of exponential functions, the ^'th eigenstate \^k)^ corresponding to the 
eigenvalue co^, is (up to an overall phase) 

V27+1 .in ((p + sffl) 



sm 



3. SU (3) phase operators 

3.1. Polar decomposition for SU (3) 

A basis for the complex extension of the Lie algebra u(3) is given by the nine operators 
{Qj : ij= 1,2,3}, with commutation relations 

[Qj ,Cki] = Sjk Qi - 5ii Ckj. (3.1) 

The complex extension of 5u(3) is obtained by restricting the u(3) operators to {Qj : / ^ j} 
and including two traceless linearly independent diagonal operators hi and h2 that determine 
a Cartan subalgebra. A convenient choice of the latter is 

/zi=Cii-C22, h2 = C22-C33. (3.2) 

If one uses the boson realization: 

Qj = d\d j , Qj = C]. , (3.3) 

the ^(A + 1)(A +2) -dimensional set of harmonic oscillator states y = {|ni^2^3),^i +^2 + 
n3 = A}, is left invariant under the action of Qj and is a basis for an irrep of 5u(3) usually 
denoted (A,0). The eigenvalues of hi and h2 are directly related to population differences 
between levels 1 and 2, and 2 and 3, respectively. 

The information in equation ( |3.1| ) is conveniently displayed using a root diagram ||47| : 
to every 5u(3) generator Cfj we associate the pair (x, of integers defined by 

[hi , Qj] = X Qj , [h2,Qj] = y Qj , (3.4) 
and the root vector a = xai + j a2, with basis ai , a2 having Cartesian components 

V2 V6^ 



ai 



The root diagram for 5u(3) is sketched in figure[T] Note that every generator is thus associated 
with a root vector and the diagonal operators are associated with vectors of length zero. 

The commutation relations are given (up to a sign) by addition of the corresponding root 
vectors. If we label the operator Qj by its root vector, Qj ia, then we have 

{e^^^p , if a + j3 is a root, 
(3.6) 
0, otherwise. 

The root diagram neatly shows, for instance, that [C23,Ci2] is proportional to C13 in 
accordance to the vectorial addition of the appropriate roots. 

Similarly, the weight diagram is a pictorial representation of the basis states of an irrep. 
The weight w = xw^ -\-yw^ of 3. basis state |w) is a vector with components related to the 
eigenvalues of the diagonal operators: 

hi\w)=x\w)^ h2\w)=y\w). (3.7) 
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C31 C32 



Figure 1. The root system for the complex extension of 5u(3), showing (in yellow) the two 
fundamental positive roots. 




= -=,-= , w2=lO,W=|, (3.8) 



The fundamental weights and have Cartesian coordinates 

1 1 

so that {w^\cCj) = 5ij. A generator Ctj associated with the root a acts on a weight state |w) by 
translation on the hexagonal grid: 

{I a + w) , if a + w is a weight, 
(3.9) 
0, otherwise. 

The operators Cfj acting on the ^ (A + 1) (A + 2) -dimensional space ^ can be represented 
by matrices. A polar decomposition of these matrices is given by 

Cij=Eijbij. (3.10) 



The operator Dij = ^ CJjQj is non-negative definite while Efj will be constructed as unitary 

matrices, with En = E22 = £33 = 1. It is easily shown that [hi^Dfj] = [h2^Dij] = and also 
that 

Dn-D2i=hu ^23-^32 = ^2, Dl^-Dj,=h^h2. (3.11) 
The realization of 5u(3) that optimally displays the polar decomposition is the 5u(3) 
analogue of (2.12). A coherent state 91, t92, 92) for the irrep (A,0) is obtained by group 
action on the highest weight state; this highest weight state is the boson state |A,0,0). A 
general quantum state |^) is then represented by a function on ^ SU {3) /U (2) 

The two azimuthal angles (pi , (p2 control the relative phase between these populations, while 
two polar angles t^i,t32 mix the number of excitations in each mode. For calculational 
convenience the latter are chosen to simplify the coherent state representation of 5u(3) 
elements, given by p9| 

hi ^r[hi) = -i- — , h2^T{h2) = -i^ — , 
C„^r«;„).ie«*---)(l + ,5^-,5^), (3.13, 
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Table 1. Notational details connecting the numbering of states, their weights, their boson 
representations, and the expression of these states in terms of the angles (p\,(p2- In the polar 
representation we omit for simplicity a factor 2k. 



State 


Weight 


Boson state 


Polar state 


|1> 


(1,0) 


1100) 




|2) 


(-1,1) 


1010) 


g/(-(Pl+(P2) 


|3) 


(0,-1) 


1001) 





&,«r,C„) = le-*.-*.(^ + .g|- + 2,g|-). 

This realization provides an obvious decomposition of the 5u(3) raising operators. Much Hke 
5u(2), these operators act in a natural way on the Hilbert space spanned by the exponential 
functions +1^292) | Again, F is not Hermitian although it is equivalent to a Hermitian 

representation, as indicated in Reference [39J or as can be shown following the procedure of 
|AppendixA 



3.2. SU (3) phase operators for the (1,0) representation 

We consider first the three-dimensional representation (1,0), spanned by the boson states 
{|100), |010), |001)}, where \n\n2n^) denotes a state with population nt in level /. This has 
been worked out in detail in f48l and from a different perspective in |49| . 

The components of a weight are related to population differences by x = ^1 — 

n2jy = n2 — n^,. Explicitly, the weights of the basis vectors are {(1,0), ( — 1,0), (0, —1)}, 
respectively. 

In this representation, we have the matrix realizations 



(3.14) 



' r 


1 


0^ 









Cl2= 










1 


i 


V 





) 















^ 









i 












• 







/ 










C23 



The rank of Cu and C23 is one, but their dimension is three, which implies that the polar 
decomposition is (again) not completely specified. Indeed, one finds that the most general 
unitary E12 and E23 consistent with the matrix realization of Cu and C23 are 

Ei2 = \ a b I , aa^'^bb'' = 1, 

(3.15) 




£23 = I 1 I , cc^^dd'' 
J* 

Here, we have already restricted Efj to be unitary so that a Hermitian phase operator can be 



properly defined. The issue is now to fix the unknown parameters a, Z?, c and J in ( 3.15 ). 
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3.2.1. SU (2) -invariant solution. The subset {Ci2,C2i,/ii} of generators spans an 5u(2) 
subalgebra of 5u(3). The boson states |100) and |010) form a two-dimensional 5u(2) 
subspace; the boson state |001) is an 5u(2) singlet. Thus, one way of fixing the unitary 
matrix £"12 is to require that E12 preserve the multiplet structure of this 5u(2) subalgebra. This 
gives a = — 1 , /? = 0, so that 

£12= -1 . (3.16) 

V 1 ; 

The phase of a is inessential and has been chosen for convenience. 

In the same way, the subset {C23,C32,/i2} spans a different 5u(2) subalgebra, and we 
may also require that £"23 act within a multiplet of this subalgebra. This, in turn, implies 




£23 = 1. (3.17) 

Vo -1 oy 

For this solution, one obtains, from Efj = q^^'J , 

/0-10\ 71 f ^ ^ ^ \ 

912 = /^ 1 0, ^3 = /^ -1 . (3.18) 
^\0 00/ ^\010/ 

The notable feature of the matrices £12 and E23 is that they do not commute. 

3.2.2. The complementary solution. An alternative choice of Eij is obtained by using a 
different line of argument. For SU{2), the phase operator is thought to be complementary 
to the population difference h. The generalization of this complementarity -based definition to 
SU (3) can also be achieved for the irrep (1,0). 

We recall the definition of the generalized Pauli matrices l|50-^|. Let 












I 


1 




(• 






























/ 



Z = /zi=|Oa)^0|, X=|0 0(0^|. (3.19) 
Then 

l^Z^ = a)^^Z^X^ (3.20) 

where G Z3 and 0) = exp(27r//3). The subset {X^Z^} of generalized Pauli matrices are 
elements of the finite Pauli subgroup of SU{3) containing 27 elements and described 



elsewhere |53 54j|. {X^Z^} also forms a basis for the 5u(3) algebra, so that we can expand 
£12 as 

En = Y.akiX''Z^. (3.21) 

k£ 

A simple analysis shows that it is indeed possible to obtain a complementary solution, in the 
sense that 

hiEi2 = co^Enh , (3.22) 

which then forces 

Ei2= \ e'P \ . (3.23) 
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A similar analysis for E23 and h2 produces the unitary solution 
/ e'> \ 

£23 = 1. (3.24) 
\ e-^^ 0/ 

If, additionally, we insist that the Ei/s commute or, equivalently, 913 = <pi2 + <P23, we 
obtain the conditions £"13 = E12E23, which implies 

j8 + 7 = 0, ±271, ... 2j3 - 7 = 0, ±2;r, ... - j3 + 27 = 0, ±2;r, . . . , (3.25) 

wherefrom we get 

3j3 =0,±2;r,... , 37=0,±2;r,... . (3.26) 

The solutions are found by fixing either j8 and deducing 7, or visa versa. The simplest 
nontrivial solution is found by choosing jS = 2;r/3 and 7 = — 2;r/3. This produces 

/010\ /Oco^OX /001\ 

En=\ (0 , £23 = 1 , £13 = 0)2 . (3.27) 
\ 0)2 00/ \a)00/ \0a)0/ 

all of which are elements of the generalized Pauli group . 



3.3. SU{3) phase operators in the (A,0) representation 



We now seek to generalize our discussion to irreps other than the simplest (1,0). The 
d = ^(A + 1)(A +2) states of the form |^i^2^3), where each nt is a non-negative integer 
subject to the condition ^1 + ^2 + ^3 = are transformed into one another under the action of 
any 5u(3) generator and thus form an irrep of dimension d. The highest weight is (A,0) and 
the highest weight state is |A00). We can again use the polar realization of equation (3.14) 
and express basis states of (A,0) in terms of exponentials. Some closely related material has 
been presented in | |55| . 

Consider for instance the case A = 2, where the dimension of the space is 6. The 
operators C12 and C23 have matrix representation and polar decomposition of the general 
form 



/0V20000\ 
' 1 ^ 







10 


0, 



■ E12 



/O 0000\ 
' V2 ^ 









10 

1 

0. 



(3.28) 



C23 


















1 





























1 











V2 















■ E23 



/ooooo 0\ 
'00000 ^ 

1 

ooooo 

.0 1 . 

Vooooov^/ 



with 



E12 ■ 



/O 1 0\ 

'000100^ 

1 

* * * 

* * * 

* * * / 



^23 



/**0*00\ 

' 1 ^ 

* * * 

1 

1 I 

\**o*oo/ 



(3.29) 



and the asterisks denoting here undetermined elements. 
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Figure 2. The construction of 5u(2)-invariant phase operators for irreps of the (A,0) type: 
equilateral triangles are turned into cones. The "rings" are made from weights in su(2)- 
invariant subspaces. 



3.3.1. SU {2) -invariant solution. One way of fixing E12 and E23 so they are unitary is to 
directly generaUze the prescription of ( 3.16 ) and ( 3.17 ) and complete the matrices in an 
5u(2)-invariant way. The 5u(2) -invariant solution, which always exists, can be obtained 
by transforming 5u(2) -invariant strings of weights parallel to a root into a circle, thus 
transforming the equilateral triangle formed by the weights into a cone, as illustrated in figure 
[2] The tip of the triangle is an 5u(2) singlet while the base is made from the longest 5u(2) 
string of weights. 



For the case of (2,0), we obtain 



E12 



/ 1 o\ 
1 » 
00010 
-100000 
01000 
\ 00001/ 



E23 



/ 1 
'001 
1 
000 
000 
\o 



0\ 
^ 

1 
1 
-10 



(3.30) 



This 5u(2) -preserving solution does not produce commuting matrices: the phases are not 
additive. This remains true for (A, 0) -type of representations, where the weight diagram is an 
equilateral triangle with A + 1 states on each side. 

One could also search for a complementary -based solution. Now, the (2,0) irrep 
of 5u(3) decomposes into a sum of two three-dimensional irreps of p^, but the resulting 
matrices are incompatible with a polar decomposition of C12. Indeed, one shows that the most 
general solution to the polar decompositions compatible with equation ( 3.29| ) cannot produce 
commuting matrices. This statement remains true for higher dimensional irreps of the type 
(A,0), provided that A is finite. 



3.3.2. Measuring non-commutativity. Quite generally the phase operators do not commute. 
To quantify the amount by which, say, £12^31 fail to commute, we introduce the matrix norm 
\\Mf= Tr (M'^M) and define 

M = Ei2E3iE-2^E^i - t. (3.31) 

Obviously, M should be the zero matrix if £12 and E31 commute. To compare values of 
||M|p for different irreps (A,0), it is convenient to normalize the length of M by dividing by 
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Ei2 




Figure 3. Left: Illustrating how the operators E12 and ^31 can commute. Red line: ^31 £'12 
acting on a state. Green line: EnE^i. Since the matrix elements of Eij is always 1, the 
operators commute in the case illustrated here. Right: Illustrating how the operators £'12 and 
£31 can fail to commute when one or the other (or both) annihilate a state. Here, the initial 
state is killed by £31 . 



the dimension ^ ( A 
solutions, we find 



1)(A + 2) of the irrep (A,0). When this is done for the 5U (2) -invariant 



\\M\\' 



•l)(A+2) 



+ !)-!] 
•l)(A+2) 



(3.32) 



This expression can be understood as follows. The action of £"12 commutes with the action of 
£31 when the action of either is non-zero, as illustrated on the left of figure [3] On the other 
hand, £12 and £31 do not commute if one or the other acts on some suitable "edge" state killed 
either by £"12 or by £31. This is illustrated on the right of figure [3] 

States killed by £12 or £"31 are always located on the edge of the weight diagram for the 
irrep (A,0). Any edge contains (A + 1) states. There are two problematic edges (one for £12 
and another for £31) so the number of problematic states is 2(A + 1). Since these two edges 
have a single state in common, the number of problematic states, adjusted for double counting, 
is just 2(A + 1) — 1. The overall multiplicative factor of 2 comes from the calculation of the 
trace, and the denominator is clearly just the normalization factor. 



3.3.3. The A ^ 00 limit and its solution. Equation ( 3.32 ) shows that, for large A, the amount 
by which £12 and £31 do not commute goes like ~ A~^: for large A, the 5u(2) -invariant 
solutions commute and the phases 912 and 931 become additive. It is also clear that, as 
A ^ 00, the edge states become progressively displaced to infinity. In this limit, the finite 
triangular lattice of the weight diagram becomes a simple two-dimensional hexagonal crystal 
lattice. 



In the polar realization ( |3.14| ), the generators acting on states having finite weight 
simplify to 

d 



h^ 



d(pi 



h2 



(3.33) 



C12 - iAe^'(2^i-^2) ^ C23 ^ ^Ae-^'(^i-2^2) ^ 



Ci3-^Ae^'(^i+^2). 
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In this limit, F is Hermitian. The rescahng Cij CijjX leads to commuting ladder operators. 

In the A ^ oo limit, the phase operators Eij act unitarily on every state of the form 
^i{n(pi^m(p2) ^j^i^ j^ j^ finite integers. The common eigenstates of £12 and £"23 are 

4. Concluding remarks 

The general prescription provided for 5u(3) can be extended to $u{n). For instance, let us 
look briefly at 5u(4). The 12 roots corresponding to {Qj : i^ 7 = 1, . . . , 4} are located at the 
vertices of cuboctahedron, which is the intersection of a cube and an octahedron. There are 
three diagonal operators, represented by three roots of length located at the center of the root 
diagram: = — Q+i^y^+i. Using again the boson realization Cij = djdj and boson states 
\n1n2n3n4), the diagonal operators correspond to population differences between consecutive 
levels. 

For such boson states, the weight diagram is a tetrahedron. Each slice parallel to a 
fundamental root of 5u(4) is an 5u(3) subspace. In particular, the action of some ladder 
operators will be undefined on one edge of the tetrahedron as the states on this edge are killed 
by them. Two polar operators will thus fail to commute when they act on states in some 
specific edge of the weight diagram. 

There are (A + 1)(A + 2)(A + 3)/6 boson states of the form \n1n2n3n4) with ^1 + ^2 + 
^3+^4 = A. There are (A + 1)(A + 2)/2 states on each edge of the weight diagram. Two 
adjacent edges intersect on a line containing A + 1 states, and they have one point in common. 



Defining M as in (3.31 ) for any pair of non-commuting roots and their phase operators, we 
found that, for 5u(4), 

^ [2x(A + l)(A + 2)-(A + l)-l] 
i(A + l)(A+2)(A + 3) ^(A + l)(A+2)(A+3) 

Our interpretation is thus that the non-commutativity of phases is an "edge" effect. Since 
the number of points on an edge of the weight diagram grows with a rate ^ A"~^ while the 
number of states in an 5u{n) representation grows like ~ A"~^ phases operators constructed 
so as to preserve the 5u{n — 1) subalgebras of 5u{n) will commute in the large A limit. 

As we have shown, in general, one cannot expect that phase operators will necessarily 
commute: they will commute in the limit of large representations where A ^ 00, a limit which 
corresponds to a contraction of 5u{n) for which F becomes Hermitian. 

We have not investigated in details the possibility of constructing commuting solutions 
which satisfy the complementary conditions for the special case of irreps of the type 
(1,0, ...,0) of 5u{n). However, it is probable that fully complementary solutions do not 
always exists. To find phase operators that are pairwise complementary is closely related to 
the existence problem for mutually unbiased bases; likely, when ^ is a prime, some elements 
of the generalized Pauli group are compatible with the polar decomposition of raising and 
lowering operators. When ^ is a power of a prime, it is not clear if the polar decomposition 
is compatible with the requirement of complementarity. When n is composite, the situation is 
even less clear as the construction of mutually unbiased bases remains an open problem. 

In conclusion, the coherent state representation used in this paper exhibits two nice 
features particularly relevant to the analysis of phase operators. First, we have a geometrical 
interpretation of them in terms of azimuthal angles related to relative phases and associated 
with the ladder action of the appropriate generators. Second, the "exponential part" of the 
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realization naturally provides the unitary part of the polar decomposition of the generators. 
Note that the similarity transformation that maps the original F realization into a Hermitian 
one simply rescales the diagonal entries of the polar part of the matrix representation of the 
generator and thus has no effect on the interpretation of phase part of the representation. 

The kind of coherent state representation used in this paper would seem to form a 
natural gateway into understanding phases in system described by algebras other than $u{n) . 
Certainly the geometrical structure of coherent states should allow one to interpret the 
parameters of such a coherent state representation. 

The authors would like to thank Laura Toppozini and Alberto G. Barriuso for numerical 
work on some aspects of this problem. This work was supported in part by NSERC of Canada, 
the Grants FIS2008-04356 and FIS201 1-26786 of the Spanish DGI and the UCM-BSCH 
program (Grant GR-920992). 



Appendix A. Making the representation F Hermitian 



In this appendix we provide some mathematical details on the representation F. The 
discussion in this section is facilitated by denoting the representation of elements of 5u(2) 
given in (2.2) by 7. In this notation, Eq.( |2.1Q[ ) takes the form 

rixm ^ [r{x)wu{(p) = {xj\RyWR,{(p)y{x)\^>) (A.i) 

The representations 7 and F are clearly isomorphic. 

Since for 2j integer every representation of 5u(2) is equivalent to a Hermitian 
representation, there must be an intertwining operator that takes the non-hermitian 
representation F in ( |2.12[ ) to the hermitian 7. 

Following reference | 56 1 we seek a similarity transformation 

r{x) = j^-^r{x)j^ (A.2) 

to bring F to its Hermitian form. Note that, under the inner product (2.13), F(/z) is already 
hermitian so commutes with F(X). Next, from the requirement that 7(^+) = 7^(^-), we 
find 

J)^'^r\e-){J(r-y = je-^T{e^)Jf (A.3) 

or J(fJf^T\e-) =F(^+)^^^ (A.4) 

Thus, the positive hermitian operator S = is the intertwining operator for which 

^F'^(^_) =F(^+)^. (A.5) 

Let F be the extension of F to the infinite dimensional space spanned by the exponential 
functions: {e^^^ : 2m G Z}. On this space we also have 

~ST\e-)=T{e^)S, (A.6) 

from which we obtain the recursion relation 

5m+i ( J + m + 1) = 5^ ( J - m) (A J) 

Because F(^_)e~^^^ = F(^+)e'^^ = 0, it is clear that, if ^{(p) is a function over the 
subspace of exponential functions spanned by {e^^^^—j <m< j}, then S^^{(p) is also in 
that subspace. Thus, we can write 5^ = 5" • 11, where 11 projects from the infinite-dimensional 
space of exponentials to the 5u(2)-invariant subspace {e^^^, —j<m< j}. The intertwining 
operator S is then the positive part of S, and hence invertible. The similarity transformation 
is the positive hermitian square root of S. is defined on {e^^^, —j<m< j}, and can 
be found recursively. 
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One may verify that the diagonal matrix with elements 

if -j<m<j, 



{j + m)\{j-m)\ 



(A.8) 



satisfies the required recursion relation. 

By simple inspection one can backcheck the formalism and verify that this "square root 
of a binomial" matrix, which we write as 



(27)! 



■ \jm){jm\, 



m=-j 



{j + m)\{j-m)\ 

satisfies [J(fJ{h)\ = [^,r(e_)r(e+)] = [^,r(e+)r(e_)] = 0, and is such that 
y{h) = J(f-^Y{h)J(f = Y{}i), 



(A.9) 



(A. 10) 



Y{e+)=Jf-'ne+)je = y{e-y. 



One immediately checks that the resulting 7 is indeed the standard Hermitian realization with 
action given in p.2|). 



Having established that T is equivalent to \22) , let us compare some aspects of their 
respective polar decomposition. The matrix realization of r(^_) is factored as 

T{e.) = T{E)nb), 



where T{b) = [r(e_)'l'r(e-)]'/2. Following equations (2.15), we find 

j 

m=-j 



(A.ll) 



(A. 12) 



r{E)= £ \j,m-l){jm\, 

m=-j 

In the same manner, using Y{e-) = 7(£) y^^/^X^^TK^-)' it established that 
r{E)= f \jm-l){jm\, 



(A.13) 



(A. 14) 



Equations (A. 14) and ( A.13 ) are identical, but neither is completely defined because the rank 
of the matrix is one less than the dimension 2j of Vj: the matrix element of E calculated 
between the lowest weight state \jj) and any other state \jm) is not determined. We can make 
E into a unitary matrix by taking m modulo 2j + 1, and thus going from the finite line to 
the circle. If we also impose the same cyclicity condition for T{E), we find T{E) and y{E) 
coincide. With the cyclic boundary conditions, the final result is: 

£ = ^ |7,m - l){jm\ , -j<m< j, mod(2j + 1) . (1.15) 

m 

A direct extension of this argument applies to the Hermitian extension of the analogous 
representation ( 3.13| ) for su(3). Details on the construction of the matrix for 5u(3) can be 
found in reference ||39J. 
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